Can liquid metal surfaces have hexatic order ? 
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We propose that extended orientational correlations can appear at the surface of supercooled heavy 
noble liquid metals, due to the same compressive forces that cause reconstruction of their crystal 
surfaces. Simulations for liquid Au show a packed surface layer structurally akin to a defected 2D 
triangular solid. Upon supercooling, the density of unbound disclinations decreases as expected in 
Nclson-Halperin's theory. It extrapolates to a hexatic transition about 350 K below melting, where 
a sharp growth of orientational correlation length and time is also found. True hexatic order is 
preempted by surface-initiated recrystallisation. 

PACS numbers: 68.10.-m, 68.35. Rh, 61.25.Mv, 82.65.Dp 



-a 
c 

o 
o 



> 

m 
O 

On 
-i— > 



i 

G 

O 

o 



X 



Heavy noble metals such as Au, Pt or Ir have a gen- 
eral tendency to reduce their lateral interatomic dis- 
tance at the surface, relative to the bulk. This tendency, 
whose microscopic origin has been traced to a relativistic 
intra-atomic s-d energy shift (!]], accounts for the zero- 
temperature reconstructions of the crystal surfaces. Re- 
markably, the lateral surface compression does not heal 
out with temperature. On the contrary, it is seen to per- 
sist, and in fact to increase, as temperature rises. For 
example, the lateral first-layer atomic density of Au(lll) 
and Pt(lll) has been found to increase, relative to the 
bulk layer density, from 1.04 to 1.07, and from 1.0 to 
1.05 respectively, between room temperature and 1000 K 
Our present rationale for this somewhat surprising 
behavior is that at high temperatures the substrate peri- 
odic potential felt by the first layer atoms is progressively 
reduced by thermal fluctuations and by a large outwards 
expansion, so that shrinkage is less hindered by epitaxy. 

In the liquid metal surface, which constitutes our 
present concern, the periodic potential is absent alto- 
gether. High lateral densities, now unhindered, can be 
expected to persist, in conjunction with strong surface 
layering, at least for temperatures which are not too high. 
Close to the melting point, X-ray data have revealed very 
pronounced layering oscillations in the density profiles 
normal to the surface in liquid metals such as Hg and 
Ga |Q, implying first layer lateral densities larger than 
the bulk average. Based on simulations, we pointed out 
some time ago || that still larger layering oscillations, 
and a higher surface lateral density, should be expected 
for the heavy noble liquid metal surfaces, boosted by the 
compressive tendencies just discussed. In addition, we 
found that supercooling makes the effect even stronger, 
with a nearly 2D-crystalline surface layer. This natu- 



rally raised the question, which we resolve in this Letter, 
whether it could be possible for hexatic, orientational 
ordering, or at least for very extended orientational cor- 
relations, to appear at the free surface of a heavy noble 
liquid metal, such as Au, Pt, or Ir, under suitable super- 
cooling conditions. 

The original suggestion made by Nelson and Halperin 
H (extending earlier work of Kosterlitz and Thouless [Q), 
and pursued by Young [|| and others || , that in 2D an in- 
termediate hexatic phase could appear between the solid 
and the liquid, led to a wide search in a variety of sys- 
tems ||. Beautiful realizations have been achieved such 
as those in 2D colloidal syste ms Jl0[ , magnetic bubble lat- 
tices [[ll| , and vortex lattices |12[. On the other hand, the 
2D simulation studies have been exceedingly controver- 
sial . An important clue which emerged more recently 
appears to be the choice of thermodynamic ensemble, 
with a strong preference of hexatic phenomena to show 
up in ensembles where density changes are naturally al- 
lowed Jig] , In grand-canonical conditions, the binding of 
disclinations taking place at the fluid-hexatic transition 
is facilitated, since it does not rely on the slow, canonical 
diffusion of defects. For instance, spontaneous creation 
of a new disclination near a pre-existing one, can lead to 
a much faster appearance of the dislocation bound state. 

In the liquid metal surface, atom exchange between 
layers can be extremely effective, with typical atomic dif- 
fusion times of order 100 ps between first and second 
"layer" for Au at the melting point |18| . This free inter- 
change of particles with the liquid substrate makes the 
outermost layer of the liquid metal a 2D system which 
is naturally "grand-canonical" on a very fast, nanosec- 
ond time scale. This is fast enough, in particular, that 
even ordinary, canonical Molecular Dynamics (MD) sim- 



1 



ulation should be able to equilibrate orientational corre- 
lations without problems. 

Moved by these motivations, we conducted very exten- 
sive MD simulations of the liquid gold surface. We used a 
slab geometry with two free surfaces, a thickness of about 
~ 60 A and a lateral size of ~ 120 A, for a total of 50000 
atoms. A second, smaller system with the same thickness 
but a lateral size of ~ 80 A was also studied in order to 
monitor size effects. Periodic boundary conditions were 
required parallel to the slab. Atoms interacted through 
the glue model hamiltonian which is well documented for 
this system, providing a good description of many prop- 
erties of bulk solid and liquid Au, of the melting tempera- 
ture, (T m — 1320 — 1350 K, against an experimental value 
of 1336), of all the main solid surfaces including their de- 
tailed energies and reconstructions, and of the liquid sur- 
face tension at the melting point Jl4||| . Constant-energy 
runs were used to generate atomic trajectories. Data 
were taken from 2000 K, well above T m , to 950 K, well 
below. Data below T m clearly refer to a supercooled liq- 
uid slab, whose lifetime we nevertheless found to exceed 
1 ns down to T - 1000 K. Below 1000 K, the lifetime 
dropped, and it became impossible to equilibrate the liq- 
uid, due to (111) recrystallization rapidly nucleating at 
both slab surfaces. Above 1000 K, where the system re- 
mained liquid, we worked at temperature intervals of 50 
K up to T m , then more sparsely. At each temperature, 
typically 30000 steps (one time-step ~ 10 -14 s) were used 
first for equilibration. Then, 20 uncorrelated configura- 
tions (40 for the smaller system) were generated from 
successive 10000-step runs, and subsequently analysed. 

The first step of the analysis consists in identifying the 
surface topmost, "first-layer" atoms in each configura- 
tion. That was done using a prescription which proved 
to work well in the past |5[. By representing atoms as 
spheres of radius R, we define surface atoms those which 
remain fully visible, i.e., totally unshadowed from above 
the slab. The optimal R value, corresponding to atoms 
in the first density profile peak in Fig. 1, was found to 
be 1.8 A, reasonably close to the atomic radius. Fig. 1 
shows a Delaunay triangulation map fl5f| for the surface 
atoms, where disclinations (five- and seven-fold) are pin- 
pointed. At high temperature there is a high concentra- 
tion of defects, both bound and free, while the density 
profile surface peak is relatively weak. Conversely, the 
peak is very strong at low temperatures, where the map 
shows large triangular domains, and a reduced density of 
free disclinations rif. Theory predicts Q a behavior of 
the form 
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as the fiuid-hexatic transition temperature Ti is ap- 
proached. We find that this equation fits very closely 
our simulation data, with a ~ 2.2 K 1 / 2 and T t ~ 974 K 
(Fig. 2). This suggests that the liquid Au surface would 
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FIG. 1. Delaunay triangulation maps for the surface 
atoms of simulated liquid Au (larger system) at T = 1450 K 
(above T m ) and T = 1000 K (supercooled). Fivefold (black) 
and sevenfold (white) disclinations are pinpointed. The two 
corresponding density profiles along the surface normal are 
also shown. The first peak correspond to first layer atoms in 
the maps. 
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FIG. 2. Density of free disclinations (of both signs) for 
the liquid metal surface as a function of temperature. Full 
diamonds: larger system; open circles: smaller system; Solid 
curve: fit to Eq. 
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undergo, in its supercooled metastable state, a genuine 
fluid-hexatic transition, if one could somehow prevent re- 
crystallization from taking place. 

Encouraged by this result, we turned to a comparison 
of positional and orientational spatial correlations, and 
of their dynamics, versus temperature. The static orien- 
tational correlation function is defined as 

ff6 (r, * = 0) = <0 6 *(r, t = 0)O J 6 (0, t = 0)) (2) 

where the orientational order parameter Oq is, at surface 
site i 

1 m 

°6 = ^E ex P( 6 ^)- ( 3 ) 

1 3=1 

The sum is restricted to the iV, Delaunay surface neigh- 
bors of atom i, and 9ij is the angle formed by the ij 
bond (projected on the xy plane) with the x axis. Posi- 
tional correlations are defined similarly, using the surface 
atom local 2D density as the order parameter. Correla- 
tion lengths £ p and ^6 are obtained by fitting the en- 
velope of the computed correlation functions to an ex- 
ponential decay, g(r) — 1 ~ exp(— r/£ p ). Fig. 3a shows 
that whereas above 1200 K positional and orientational 
surface correlation lengths are close, and about 5 A, be- 
low this temperature the latter take off while the for- 
mer do not. A fit of £e(T) to theory, which predicts 
£ 6 (T) ~ exp [a/(T — X^) 1 / 2 ] yields reasonable agreement 
for our larger system, with a ~ 15.5 K 1 / 2 and Tj ~ 942 K. 
There is however an alarmingly large discrepancy with 
the parameter values of the disclination fit. The disagree- 
ment turns out to be instructive, and is explained by a 
size effect. By comparing our small and large systems, we 
observe that whereas disclination densities hardly depend 
on size, the growth of long-ranged orientational correla- 
tions is enormously hampered by size. As Fig. 3a shows, 
a decrease of lateral size from 120 to 80 A is sufficient to 
damage the rise of reducing it (at 1000 K) from 12.5 
to 7.5 A. This is, we surmise, an independent signature 
of the Coulomb nature of effective interactions between 
disclinations in the liquid metal surface. Size effects in 
two-body correlations are known to be extremely large in 
related systems such as the XY model Jig] . Fortunately 
(a proper size scaling being presently out of the ques- 
tion), the one-body disclination density is instead well 
converged (Fig. 1). We conclude that, even though & is 
still far from convergence even in the large system, the 
disclination fit parameters should safely represent the in- 
finite system. 

As theory suggests [jl7|, and as the colloidal ex- 
periments have demonstrated [ {To) , the temperature- 
dependent orientational relaxation times can provide an 
even stronger dynamical signature of hexatic phenomena. 
We computed the positional and orientational correlation 
times t p and tq by an exponential time fit of the envelope 
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FIG. 3. (a) Orientational (full diamonds) and positional 
(empty diamonds) correlation lengths. Open circles: orien- 
tational length for the smaller system. The full line is a fit 
to theory (see text), (b) Orientational (full diamonds) and 
positional (empty diamonds) correlation times. Dotted lines 
are guides to the eyes. 

of the time-dependent 2D structure factor F(K, t), where 
K is the frustrated reciprocal lattice spacing, and of the 
time-dependent orientational correlation function (r = 
0, t). Fig. 3b shows the results for the larger system. Be- 
low about 1200 K tq is found to rise dramatically (note 
the log scale) while t p does not. There is a clear orien- 
tational slowing down, which complements very well the 
static results, and also suggests that a likely experimental 
signature to be sought at the liquid metal surface might 
be dynamical rather than static. 

In spite of these growing correlations, true hexatic 
long-range order appears, in our particular system, to 
be narrowly preempted by recrystallization. Already at 
1000 K, after a lifetime of about 0.5 ns the liquid slab sud- 
denly starts recrystallizing. At 950 K, the liquid lifetime 
(starting with the 1000 K liquid configuration) drops to 
less than 0.1 ns. Recrystallization proceeds with two 
fronts starting simultaneously at the two surfaces. At 
1000K, we analysed first the metastable liquid surface, 
and then the moving recrystallization front, monitoring 
values of £ 6 an d £ p in the successive layers as a function 
of depth jl8|. As shown in fig. 4 the front is split, an ori- 
entational ordering front preceding a positional ordering, 
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FIG. 4. Orientational (diamonds) and positional (trian- 
gles) surface depth profiles. Bottom: liquid surface, T=1000 
K. Top: during recrystallization. Note that the orientational 
front precedes the positional front. 

crystalline front. The intermediate region, where orien- 
tational correlations are extended, but positional corre- 
lations are short, is strongly layered (not shown), and 
extends for a thickness of about three layers. In fact lay- 
ers just below the surface already possess, in the purely 
liquid slab, long orientational correlations (fig. 4). This 
kind of surface wetting by an orientationally ordered film 
coexists, as we verified directly, with rough, liquid-like 
height- height static correlations (not shown), and seems 
an exclusive characteristics of a strongly layered liquid 
surface. No such orientationally ordered film is found, 
for example, in a Lennard-Jones system, either at the 
free liquid surface, or at the crystal-liquid interface p9[ . 
At the liquid metal surface, the presence of the orienta- 
tionally ordered film seems to provide the required germ 
for propagation of the crystal phase into the bulk. 

In conclusion, the supercooled surface of liquid Au, 
and most likely also that of other heavy noble metals 
such as Pt or Ir, is predicted to develop an incipient 
long-ranged orientational order described by the Nelson- 
Halperin theory. The order extends a few layers below 
the surface. Transition to a well defined hexatic phase is 
preempted by a fast and effective surface-driven recrys- 
tallisation, where the orientationally ordered film pre- 
cedes the crystalline front. Our results should motivate 
future experiments, either structural, to detect the ori- 
entationally quasi-ordered film, or dynamical, addressing 
the orientational slowing down, or metallurgical, to inves- 
tigate the nature and depinning of the recrystallization 
front in these simple but exciting systems. 
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